NASA CR-1 68266 
GTD 83-10 


r " 1 



INVESTIGATION OF A 


PULSED ELECTROTHERMAL THRUSTER 


by R. L. Burton, S. A. Goldstein, B. K. Hllko, 
D. A. Tidman, and N. K. Wlnsor 


i 

GT-DEVICES, INC. 


prepared for | 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


1 

4 


NASA Lewis Research Center 
Contract NAS 3-23779 




■ 64 - 1222 ? 


03/20 


Unclaa 

42578 


4r v 






0 


1. Raport No. 

NASA CR- 168266 


2. Government Accession No. 


4. This and Subtitle 

Investigation of a Pulsed Electrothermal Thruster 


7. Author(s) 

R. L. Burton, S. A. Goldstein, B. K. Hilko, 0. A. Tidman, 
and N. K. Winsor 


9. fVforming Organization Name and Address 

GT-Devices, Inc. 

5705 General Washington Drive 
Alexandria, VA 22312 


12. Sponsoring Agency Name and Address 
National Aeronautics and Space Administration 
Washington, DC 20546 


3. Recipient's Catalog No. 


5. Report Date 

October 1 . 1983 


C Performing Organization Code 


8. Performing Organization Report No. 

GTD Report No. 83-10 


10. Work Unit No. 


11. Contract or Grant No. 
NAS 3-23779 


13. Type of Report and Period Covered 
Contract Report 


14. Sponsoring Agency Code 


IS. Supplementary Notes 

Project Manager, Shih-YIng Wang, Space Propulsion Technology Division, 
NASA Lewis Research Center 


16. Abstract 

Exhaust velocity and thrust measurements are. performed on a Pulsed Electrothermal (PET) 
Thruster using polyethylene and Teflon propellants. The results verify theoretical predictions 
of equilibrium flow in the nozzle, resulting In substantial recovery of the energy of 
dissociation and ionization. The thruster is tested in an unsteady mode (15 nsec current pulse 
and 15 cm discharge length) and in a quasi-steady mode (48 nsec current pulse and 5 cm discharge 
length). All tests are run at 2 kJ. The exhaust velocity of the propellant mass exiting during 
the current pulse is measured with two types of time-of-flight probes, and the impulse bit is 
measured on a thrust stand. It is inferred from both theory and experiment that an additional 
amount of mass is exhausted after the pulse. The measured thrust to power ratio for polyethylene 
Is T/P = 0.10 N/kW at 21 km/sec in the unsteady mode, and T/P = .053 N/kW at 27 km/sec in the 
quasi-steady mode, where the velocities are measured by the time-of-flight probes. For Teflon 
propellant, T/P = .20 N/kW at 15 km/sec (unsteady mode) and 0.090 N/kW at 20 km/sec (quasi-steady 
mode). The discharge pressure and temperature predicted by a computational model for polyethylene 
are consistent with the measured thrust and discharge resistance. The total measured propellant 
ablation Is in excellent agreement with that given by the model, which predicts that about 2/3 
of the mass is ablated during the current pulse and about 1/3 at later times. The computed sound 
speed for polyethylene Is 9 km/sec. Indicating that the exhaust plasma (27 km/sec) is highly 
supersonic. 


17. Kay Wo rds (Suggsstad by Author(s) ) 

Electric Propulsion 
Electrotnermal Thruster 
Thruster Performance 


18. Distribution Statamant 

Unclassified-Unlimited 


19. Security Qmlf, (of this rtport) 

20. Security CJascif. (of this psgt) 

21. No. of Pagn 

22. PHct* 

Unclassified 

Unclassified 

40 



* For sale by the National Technical Information Service, Springfield, Virginia 22161 


NASA-C-168 (Rav. 10-75) 


1 





TABLE OF CONTENTS 


Pa^e 

Title Page 1 

Table of Contents 11 

Summary Ill 

I. INTRODUCTION 1 

II. DESCRIPTION OF PET THRUSTER OPERATION 3 

III. DESCRIPTION OF EXPERIMENTAL APPARATUS 10 

IV. EXPERIMENTAL RESULTS 25 

V. DISCUSSION OF RESULTS 28 

VI. CONCLUSION 39 

VII. APPENDIX I-REPORT DISTRIBUTION LIST .40 

VI I I. REFERENCES 43 





SUMMARY 


The PET Thruster is an electric propulsion device which 
consists of a capillary-shaped pressure chamber with electrodes 
at each end. The cathode end is closed, and the exhaust 
electrode (anode) incorporates a supersonic nozzle. The pressure 
chamber inner wall is lined with an easily ablated, low molecular 
weight material such as polyethylene or Teflon. The pressure 
chamber is a high strength dielectric, capable of withstanding 
plasma pressure pulses up to 104 atm. The high pressure develops 
a jet of high velocity plasma from the nozzle, which provides 
rocket thrust. 

The electrical energy store consists of a high voltage 
capacitor or pulse forming network, capable of storing a few 
kilojoules at several kilovolts. The impedance of this source is 
matched to the discharge resistance R, for efficient coupling. A 
switch with resistance R S <<R permits over 90% of the stored 
energy to be transferred to the discharge in a single pulse. The 
width of the current pulse Is typically 10-100 fisec, with a peak 
current In the 20-100 kA range. Terminal voltage is 3-10 kV, giving a peak 
power dissipation on the order of 100 MW. 
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During a current pulse, 1 2 r heating in the discharge chamber 
creates a plasma which radiates as a black body to the wall. The 
wall material is evaporated and is heated to a temperature of a 
few eV. The chamber pressure is typically ~ 10, 000 atmospheres, 
so that the plasma Is highly collisional and optically thick. 
For steady conditions in the chamber, the pressure drives a 
steady plasma flow through the supersonic nozzle. Because of the 
high gasdynamlc pressure this nozzle flow is predicted to be in 
equilibrium. Thus, the efficiency of the PET Thruster is 
expected to be high. 

Exhaust velocity and thrust measurements described in this 
report are performed using poiyethylene and Teflon propellants. 
The results verify the theoretical predictions of equilibrium 
flow in the nozzle, resulting in substantial recovery of the 
energy of dissociation and ionization. The thruster is tested in 
an unsteady mode (15 n sec current pulse and 15 cm discharge 
length) and In a quasi-steady mode (43 //sec current pulse and 5 
cm discharge length). All tests are run at 2 kJ. The exhaust 
velocity of the propellant mass exiting during the current pulse 
is measured with two types of tlme-of-f 1 i ght probes, and the 
Impulse bit Is measured on a thrust stand. It is inferred from 
both theory and experiment that an additional amount of mass is 
exhausted after the pulse. The measured thrust to power ratio 
for polyethylene Is T/P * 0.10 N/kW at 21 km/sec In the unsteady 
mode, and T/P » 0.53 N/kW at 27 km/sec in the quasi-steady mode, 
where the velocities are measured by the ti me-of-f 1 1 ght probes. 



For Teflon propellant, T/P = .20 N/kW at 15 km/sec (unsteady 
mode) and 0.090 N/kW at 20 km/sec (quasi-steady mode). 

The discharge pressure and temperature are predicted by a 
computational model for polyethylene, and are consistent with the 
measured thrust and discharge resistance. The total measured 
propellant ablation Is In excellent agreement with that given by 
the model, which predicts that about 2/3 of the mass is ablated 
during the current pulse and about 1/3 at later times. The 
computed sound speed for polyethylene Is 9 km/sec, Indicating 
that the exhaust plasma (27 km/sec) is highly supersonic. 
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I. INTRODUCTION 


After more than two decades of electric propulsion develop- 
ment, a requirement still remains for a system which has high 
efficiency In the 1000-2500 second specific Impulse range. The 
problem is Illustrated In Figure 1, which shows overall effici- 
ency as a function of specific Impulse for several electric 
propulsion systems.! Below 750 seconds and above 2500 seconds 
efficiencies exceeding 70% can be obtained, but In the Intermed- 
iate range the efficiency of available systems falls as low as 
30%. 


The reason for the poor performance of these systems at 
1000-2500 seconds Is primarily the large Internal energy of the 
exhaust flow, which Is highly dissociated or Ionized, represent- 
ing a "frozen flow" energy loss which Is not recoverable. Because 
the flow Is frozen by the low number of recombination collisions, 
raising the flow density several orders of magnitude raises the 
collision rate to where most of the plasma Internal energy can be 
recovered. In this report we present the results of the first 
thrust stand and exhaust flow tests of the PET (Pulsed Electro- 
thermal) Thruster, In order to check our prediction! of 
equilibrium flow In the nozzle. If true, equilibrium flow would 
permit efficient operation of the PET Thruster In the range of 
1000-3000 seconds specific Impulse. 



Figure 1 
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II. DESCRIPTION OF PET THRUSTER OPERATION 


The PET Thruster (Figure 2) is an electric propulsion device 
which consists of a capillary-shaped pressure chamber with 
electrodes at each end. The cathode end Is closed, and the 
exhaust electrode (anode) incorporates a supersonic nozzle. The 
pressure chamber Inner wall Is lined with an easily ablated, low 
molecular weight material such as polyethylene or Teflon. The 
pressure chamber Is a high strength dielectric, capable of 
withstanding plasma pressure pulses up to 104 atm. The high 
pressure develops a jet of high velocity plasma from the nozzle, 
which provides rocket thrust. 

The design of the thruster permits the propellant Inner 
diameter and length to be changed easily for testing. The 
propellant can also be easily removed for weighing or replace- 
ment. The nozzle design places the anode downstream from the 
throat In a 1 ow-temperature region, which virtually eliminates 
erosion on this electrode. 

The electrical energy store (Figure 2) consists of a high 
voltage capacitor or pulse forming network, capable of storing a 
few kilojoules at several kilovolts. The Impedance of this 
source Is matched to the discharge resistance R, for efficient 


Figure 2. PET THRUSTER SCHEMATIC 
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coupling. The switch (Figure 2) Is an Ignltron with resistance 
R $< <R , so that over 90% of the stored energy Is transferred to 
the discharge In a single pulse. The width of the current pulse 
Is typically 10-100 #isec, with a peak In the 20-100 kA range. 
Terminal voltage Is 3-10 kV, giving a peak power dissipation on 
the order of 100 MW. 

It Is possible to operate the PET Thruster In two modes, and 
results from both modes are presented In chls report. The mode Is 
determined by comparing the length of the current pulse, tp, to 
the two-way longitudinal acoustic time In the capillary, Zi/c. 
In the short pulse or unsteady mode, t p c/2*< 1, and steady flow 
Is not achieved through the nozzle. In the long pulse or quasi- 
steady mode, tp C /2*>> 1, and steady flow Is achieved during the 
high current portion of the pulse. The 2 kJ tests described 
below result in higher specific Impulse and efficiency for a PET 
Thruster operated In the quasi-steady mode. The two current 
pulses tested are shown In Figure 3. 

Physics of Operation 

The physics of the PET Thruster are most easily described 
for constant current operation.! Electrothermal heating In the 
discharge chamber at a rate I?R creates a plasma which radiates 
as a black body to the wall. The wall material is evaporated and 
Is heated to a temperature of a few eV. The chamber pressure Is 
typically several thousand atmospheres, so that the plasma Is 
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highly collislonal and optically thick. For steady conditions in 
the chamber, the pressure drives a steady plasma flow through the 
supersonic nozzle. Because of the high gasdynamlc pressure this 
nozzle flow is In equilibrium, as discussed below. Thus, the 
efficiency of the PET Thruster is expected to be high. 

Nozzl e Processes 

Without an equilibrium flow nozzle, a highly ionized plasma 
exhaust would result in high frozen flow losses. At high 
densities, however, it appears that the ionization energy can be 
recovered. The three-body recombination rate is^ 


v 3 8 8.75 x 10 ” 27 T e y 4 ' 5 Z 3 n 2 (cm" 3 ) [sec" 1 ] 

We can compare the collislonal recombination time 
based on expected plasma conditions, to the nozzle flow time 
*/u e . The smallest value of **3 occurs at the nczzle exit plane. 

or estimated values TeV 85 0-5 eV , ne s 2.5 x 1 0 1 8 ^ Ug * 2.5 x 
106 cm/sec, and nozzle length *- * 10 cm, v 3 Jt/u g as 2 x 10 7 » so 

that recombination collisions are frequent in the nozzle. The 
flow expansion process therefore remains close to Saha equilib- 
rium, and the ionization can be largely recovered. 
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Estimated Thrust Efficiency 
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We have tabulated the chemical bond energies for poly- 
ethylene, 3 in order to estimate the thrust efficiency for equil- 
ibrium flow. The relevant energies for the production of six 
Ions (2 carbon and 4 hydrogen) are summarized below: 


evaporate C 2 H 4 from wall: 0.5 eV 

break CC double bond 4.5 eV 

break 4 CH bonds at 3.5 eV 14.0 eV 

total for 6 Ions 19.0 eV 

evaporation and dissociation energy/Ion 3.2 eV 

mean Ionization energy/Ion 12.8 eV 

total 16.0 eV/lon 


We now use the Saha equation to estimate the degree of 
Ionization at the exit plane. At low temperatures the Saha 
equation can be simplified to 


a * 3.11 x 10^ p"* T fi y^e“ e 1^^eV 
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where the pressure is in MKS units. For p e = 5 x 10 4 N/m2, 
T e v * 0.5eV , e diss = 4.5 eV and e i * 12.8 eV, the Saha equation 
predicts ^Iss = .65 and <Mon 55 1.6 x 10" 4 * We can thus estimate 
the mean energy loss for the six ions as 

loss from evaporation 0.5 eV/ion 

loss from dissociation (.65 x 18.5/6) 2.0 

loss from ionization (.00016 x 12.8) 0 

loss fron thermal (3/2 TeV^ -lii 

E loss s 3.3 eV/ion 

thruster efficiency for this example Is then calculated 

n >U I 

" T *“e * eE loss /m 

The mean Ion molecular weight for polyethylene Is 4.67. As an 
example, for Eloss * 3.3 eV and ue/9 = 2500 seconds, n T = *®2. 


The 
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III. DESCRIPTION OF EXPERIMENTAL APPARATUS 


The performance of the PET thruster Is measured on a simple 
linear thrust stand (Figures 4 and 5). The measurement Is per- 
formed oy freely suspending the thruster, firing It, and measur- 
ing the recoil velocity V. Knowing the recoil mass M, th? 

Impulse bit Is measured by 


/ 


Tdt [N-s] = MV 


In order to Insure that the thruster Is freely suspended, 
the design of the thrust stand Incorporates the following ele- 
ments (see Figure 6). 

1. The thruster mass (2.5 - 2.7 kg) Is supported by a 
linear bearing which rolls on a horizontal hardened steel shaft. 
The center of gravity of the recoil mass lies directly on the 
thrust axis and In the vertical plane directly below the center 
of the bearing. Firing the thruster therefore produces simple 
linear motion with no binding of the bearing. 

2. Current Is fed Into and out of the thruster through 
flexible strips of copper window screen (Figure 6). The strips 
are pie shaped with minimal gaps between the segments to minimize 
the fringing magnetic fields which could affect the diagnostics. 
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3. The thruster Is sealed to the vacuum tank by a thin (0.4 
mm) flexible diaphragm of rubber fabric, located at the nozzle 
throat (Figure 6). This arrangement keeps the electrical connec- 
tions outside the vacuum tank, greatly simplifying the design of 
the high voltage Insulation. The thruster and diaphragm must 
sustain a load caused by atmospheric pressure, and this load is 
counterbalanced by coll springs (Figure 6). The natural 
frequency of the thruster mass-spring system Is 16 Hz, and the 
damping of the system Is low, as Is shown In the upper plot In 
Figure 7. 

The recoil velocity Is measured by a non-contact tech- 
nique, using an Inductive RF trandsucer4 coupled to an AISI 4142 
steel target. The sensitivity of the transducer Is 8V/mm, so 
that micron-level motions can be detected. The transducer output 
signal Is recorded on a digital oscilloscope (Smartscope Model 
3281) at 2 //sec/polnt and Is graphically differentiated to give 
V, the Initial recoil velocity (Figure 8). Alternatively, the 
recoil velocity can be determined by fitting the output signal to 

a waveform x(t) * x 0 e -ot s1no>t, where <*, x 0 and » can be easily 
measured (Figure 7). Both techniques give equivalent results. 

Delivered Energy 

PET Thruster performance Is based on the energy delivered to 
the discharge, determined by integration of the current and 
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voltage at the thruster terminals. Thruster current I Is 
measured with a Rogowski loop, passively integrated with an RC = 
2.8 msec integrator. Terminal voltage V is measured with a 935:1 
resistive divider. Both I and V are recorded at 2/isec/point, 
enabling the product IV to be internally computed, integrated, 
and displayed on the scopeface (Figure 9). The voltage is V = IR 
+ L T i, where Lt is the inductance of the thruster. Thus, 

fmt = fl* Rdt + i L t I 2 )dt 

Denoting Ei as the integral of the first half cycle, then 
E 1 = fmt = /I 2 Rdt = E D [J] 

where E D i s the delivered energy. Because of the small external 
resistance, Eg 1S typically 92% of the stored energy. 

Ablated Mass 


The propellant ablated mass m is determined by a weighing 
procedure. The propellant is weighed on a sensitive scale 
(Voland Model 220R), and is weighed again after firing one or 
more shots. The cathode is weighed similarly. The nozzle is not 
weighed; however, visual inspection after 100 shots on the 
nozzle has shown no visible erosion. 



Figure 9 
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Exhaust Velocity 


A number of techniques have been employed to measure exhaust 
velocity. These include a) an ion saturation probe, b) a piezo- 
electric stagnation pressure probe, and c) photodiode sensors, 
which record time of arrival of high velocity plasma at a stag- 
nation point. The velocity measured by each of these techniques 
was generally In agreement. A description of each technique Is 
presented below. 

The Ion saturation probe consists of 4 flush wires of 0.1 mm 
diameter mounted In a 90° epoxy cone. The wires are biased to 
-24V with an automotive lead-acid battery, referenced to the 
center electrode. When plasma arrives at the probe tip, a satu- 
ration current flows through the probe given by the saturation 
current density j + = 1 /4 n + e v -1 ^ flowing to the exposed area of 
the wires. Unfortunately this technique Is unreliable, because 
the flush mounted wires become covered with material after a 
small number of shots, producing loss of signal. This covering 
Is probably melted epoxy material from the probe tip. 

The stagnation pressure probe (Figure 10) consists of a 
sensitive piezoelectric transducer (PCB J.02A12, 44 p$1/volt), 
mounted on the vacuum tank centerline on a moveable support. The 
transducer tip Is covered with mylar tape for protection from the 
discharge. Typical probe signals are shown In Figure 11 for two 


19 






Figure 11 

PRESSURE PROBE RESPONSE vs POSITION 


ORIGINAL ?AGE 53 
OF POOR QUALITY 


PROBE @ 
98cm 
(psi) 


PROBE© 
54 cm 
(psi) 


THRUSTER 

CURRENT 

(kA) 


POLYETHYLENE 
5.5cm CAPILLARY 
18 kJ DELIVERED 



TIME (yj. sec) 


21 



probe positions, showing that a ti me- of-arri val can be clearly 
measured for each shot. These times are then plotted, and the 
mean slope is measured graphically to determine the exhaust 
vel ocity . 

A third time-of-arri val technique employs simple red-sensi- 
tive diodes to record the plasma flash on the tip of a 
transparent plastic stinger (Figure 10). The probe is located 
downstream from the nozzle exit plane and records a spike of 
radiation when fast moving plasma strikes the probe (Figure 12). 
Again, the arrival times for different ax*al probe positions are 
plotted to determine the exhaust velocity. 

Because the total mass of the ambient gas in the tank at 

p amb~*003 torr is a small fraction of the exhaust mass 4m, both 
the pressure transducer and the photodiode probes can be expected 
to respond to the exhaust mass without interference of the back- 
ground gas. The data from the photodiode probe predicts 
velocities consistenc with the data from the pressure transducer. 
Typical arrival time data, and the graphical measurement of 
exhaust velocity, is shown for both probes in Figure 13. 
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Figure 12 
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Figure 13 
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IV. EXPERIMENTAL RESULTS 


Two separate sets of tests were run on the PET Thruster, 
using two different current pulse waveforms (Figure 3). The 
first set of tests was run In the unsteady (short pulse) mode 
for which the thruster was driven by two 29 /if d capacitors. 
Because of the high plasma resistance (0.2Q), approximately 92% 
of the stored energy Is transferred to the discharge chamber. 

The second set of tests was run In the quasi-steady mode, 
for which the PET Thruster Is driven by a broad current pulse, of 
48 ftiec duration, (Figure 3). The pulse was again switched by an 
Ignltron. The propellant chamber used In the quasi-steady mode 
Is 5.5 cm long with a minimum resistance of .1600, so that 90% of 
the stored energy Is transferred. 

The test plan calls for a ten-shot series with each propel- 
lant, polyethylene and teflon, so that 20 shots were run In the 
unsteady mode. An additional shorter series of shots was run In 
the quasi-steady mode. Yhe data recorded, and the estimated 
precision of the measurements are: 


25 


(1) recoil motion x(t). Impulse bit Is measured to 5%. 

(2) Exhaust probe tlme-of-f 1 Ight, accurate to 10%. 

(3) I, accurate to 3%. 

(4) V, accurate to 3%. 

(5) Mass loss 4m, accurate to 6%. 

The basic data are shown below In Table III and Table IV as 
averages for each series. These results are analysed in the next 
section . 



Table III. PET THRUSTER PERFORMANCE - UNSTEADY MODE, 

15 CM CAPILLARY 



Polyethylene 

Teflon 

Delivered energy. Eg 

1.80 kJ 

1.87 kJ 

Maximum current, I M 

23 kA 

22 kA 

Pulse width, half maximum 

15 jisec 

15 n sec 

Thrust Impulse bit, /Tdt 

.20 N-s 

.35 N-s 

Exhaust velocity 



Pressure transducer 

21 km/s 

15 km/s 

Photodiode 

20 km/s 

14 km/s 

Ablated mass, capillary only 

20 mg/shot 

117 mg/shot 


Table IV. PET THRUSTER 

PERFORMANCE-QUASI- 

STEADY MODE, 

5.5 CM CAPILLARY 



Polyethylene 

Teflon 

Delivered Energy, Ed 

1.78 kJ 

1.89 kJ 

Maximum current, I M 

17 kA 

12 kA 

Pulse width, half maximum 

48 |isec 

48 //sec 

Thrust Impulse bit, /Tdt 

.095 N-s 

.17 N-s 

Exhaust velocity 



Pressure transducer 

28 km/s 

20 km/s 

Photodl ode 

26 km/s 

20 km/s 

Ablated mass, capillary only 

6.5 mg/shot 

27 mg/shot 
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V. DISCUSSION OF RESULTS 


At our present level of understanding and interpretation of 
the data, the results presented above in Tables III and IV should 
be interpreted cautiously. The thrust stand response, while 
highly repeatable, has not as yet been checked with an 
independent calibration procedure. Further uncertainty exists 
in the interpretation of the exhaust velocity and ablated mass 
measurements, as discussed below. 

The thrust stand was designed as an oscillating spring-mass 
system with low damping, and this goal has been achieved, as 
shown in the slowly decaying sinusoidal thrust stand response in 
Fig. 7. Confidence In the accuracy of the thrust stand results 
Is based on this low damping, and on the repeatability of the 
measurements of slope at t * 0. The raw Impulse bit data for the 
2 kj tests presented In Tables III and IV is shown In Table V 
below. These results show that the thruster Is operating 
repeatably, and that the measurement of the thruster recoil 
velocity does not Introduce appreciable uncertainty Into the 
final result. 
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Table V. PET THRUSTER IMPULSE BIT DATA, 2 kJ SHOTS 


A. Polyethylene, 15 ysec pulse B. Polyethylene, 48 ysec pulse 


/Tdt 


.184 N-s 

.184 

.190 

.198 

.205 

.199 

.204 

.195 

.194 

.193 


Mean .195 ±.007 N-s 


C. Teflon, 15 ysec pulse 


/ Tdt * .373 N-s 
.340 
.315 
.406 
.368 
.350 
.334 
.400 
.335 
.311 


Mean .353 + .033 N-s 


/Tdt = 


.094 N-s 

.0965 

.0966 

.0940 

.0946 

.0967 

.0958 


Mean .0954 ±.0014 N-s 


D. Teflon, 48 ysec pulse 


/Tdt « .169 N-s 
.170 
.169 
.159 


Mean .167 + .005 N-s 


P 
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The measurement of exhaust velocity gives results which are 
consistent with our understanding of the physics of the discharge. 
In the unsteady mode case (Table III), the current impulsively 
heats the discharge plasma which is then accelerated as the 
temperature decays. The measured exhaust velocities are 20-21 
km/sec for polyethylene and 14-15 km/sec for Teflon. The quasi- 
steady case (Table IV) achieves nearly steady state conditions, 
and hence a higher mean sound speed for the duration of the 
pulse. The exhaust velocities are correspondingly higher: 26-28 

km/sec for polyethylene and 20 km/sec for Teflon. The ratio of 
exhaust velocity for polyethylene Is about 1.4 that of Teflon, 
primarily due to Teflon's higher molecular weight. 

What evidence do we have that the exhaust velocities measur- 
ed correspond to electrothermally heated and accelerated propel- 
lant mass, and not to some relatively low mass precursor or shock 
wave? First is the time response of the pressure and photodiode 
probes (Figs. 11 and 12), which display a pulse width about equal 
to that of the current pulse. This suggests that the probe is 
responding to the momentum flux generated In the nozzle. Second 
Is the nearly comparable velocity measured by the pressure and 
photodiode probes (Fig. 13); the latter responding to the 
radiation emitted by exhaust mass stagnating on a blunt body In 
the exhaust stream. Third, we have a theoretical prediction of 
the exhaust velocity, described In the next section, which gives 
reasonably good agreement with the measured values. 
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Computati onal M odel 


The predicted performance of a PET Thruster capillary with 
polyethylene can be modeled numerically.! Inputs to the GT- 
Devices code include polyethylene chemistry, 5 plasma resistivity, 
capillary geometry, and the experimentally measured current and 
voltage waveforms. Because the plasma resistance is known, this 
can be used to calibrate the resistivity model used in the code, 
with the result that the plasma resistivity is close to that 
given by the Spitzer model. 6 Predicted performance curves for a 
polyethylene 2 kJ shot in the quasi-steady mode are shown in 
Figure 14, which shows the capillary pressure, temperature and 
sound speed. The sound speed c 0 the capillary is 8.5 km/sec 
during the current pulse, during which timey«1.2. Assuming 
adiabatic flow, with constant y* 1.2 in the nozzle, the exit 
Mach number is then M~4.7, corresponding to the nozzle area 
ratio of 95. At this Mach number the temperature drops by a 
f actor : ' 

^ 1 + M2 = 3.2 


The sound speed thus drops by a factor V3.2 = 1.8, to c 5 4.8 

km/s. The exhaust velocity is then u e = Me * 22 km/sec, approximately 
the measured exhaust velocity. 

More evidence for high Mach number flow lies in the pressure 
probe response (Fig. 11), which shows a maximum signal of about 2 
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atmospheres. The thruster stagnation pressure can be estimated 
from the measured Impulse bit, throat area, and pulse length? as 


P ~ -710 atm 

C F fl *'p 

for Cp~2 and tp^50 fisec. This pressure Is 1 .1 agreement with 
that predicted by the computational model (Fig. 14). If we assume 
that a normal shock sits on the pressure probe tip, then In the 
limit m 2 >> (y-i)/(y+i) the probe response p 0 y t0 the shocked gas 
Is given by? 

p o 

The probe peak amplitude of 2 atm corresponds to p oy * 350 . For 
y- 1.2, this pressure ratio Indicates M~6.3, corresponding to a 
flow velocity of 24 km/sec for c 0 « 8,5 km/sec. This value Is 
consistent with the exhaust probe tlme-of-fl Ight measurements. 

Based on the numerical model, the measured exhaust 
velocities are consistent with a highly supersonic, adiabatic 
nozzle. Consistent with this model, the pulse width detected by 
the pressure transducer (Fig. 11) remains roughly constant with 
distance from the nozzle. This non-dlspersl ve feature Is charac- 
teristic of a high velocity puff of gas with a low sound speed. 
Furthermore, the measured velocity Is consistent with adiabatic 
expansion of a polyethylene plasma with a stagnation temperature 
of 2.5 eV and a stagnation pressure of 10 3 atmospheres. The 
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stagnation enthalpy Is 340 kJ/gram for polyethylene plasma, 5 
which would then fully expand to 26 km/sec. This value Is close 
to those measured by the velocity probes. 


Thrust to Power Ratio 

ratio [N/kW] can be evaluated in the 

/Tdt _ /Tdt 
‘ /Pdt'/TOOO " E^[kJ] 


The thrust to power 

form 

} [N/kW] 


The values of T/P for the 2 kJ shots (Tables III and IV) are 
shown below In Table VI. The velocities shown are those measured 
by the velocity probes for the mass exhausted during the current 
pulse. 


Table VI. THRUST TO POWER RATIO FOR 2 kJ PULSES 



Polyethyl ene 

Tef 1 on 

Unsteady, 

.10 N/kW 

.20 N/kW 

15 n sec pulse 

021 km/sec 

015 km/sec 

Quasl-iteady , 

.053 N/kW 

.090 N/kW 

48 usee pulse 

027 km/sec 

020 km/sec 
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OF POOR QUALITY 


For the quasi-steady case, these T/P values can be used to infer 
a thruster efficiency, as discussed below. 

Thruster Efficiency 

The thruster efficiency is conventionally defined as 

_ exhaust kinetic energy _ * mu e 
n T energy delivered ~ Ep 

For a pulsed device, the exhaust kinetic energy Is 

i mu| = fi mu| dt 

where the Integral Is continued past the end of the current pulse 
length t p to account for the discharge afterglow, during whi~h 
propellant evaporation decays to zero. If the PET Thruster is 
operating in a purely steady mode, the efficiency can be found 
from u e and T /p by 


n T = 



Using this efficiency and the measured u e the energy loss per 

particle can be estimated from the equation for n T S^ven in 
section II above. 

Thus the quasi-steady data in Table VI would predict: 


35 


Polyethylene, 27 km/sec 


n T = .72 Eloss “ 6 * 9eV 

Teflon, 20 km/sec n T = .90 Eloss = 3.7eV 


While these numbers are certainly encouraging, they tend to 
overstate the efficiency, because the mass exhausted after the 
pulse at low velocity Is not included in the calculation. As a 
check on the post-pulse mass ablation and on the quasi -steady 
nature of the discharge, we refer to the ablated mass 
measurements in Table IV. These mass measurements can be used to 
define a time-averaged exhaust velocity u e b y 


The values of u e and ue for the quasi-steady data are shown in 
Table VII. 


Table VII. MEASURED u e AND PULSE-AVERAGED u e » 
QUASI-STEADY MODE 

Polyethylene Teflon 

Measured ue 27 km/sec 20 km/sec 

Mean u e = JLjTdt 15 km/sec 6 km/sec 

Am 
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The values in Table VII indicate that the thruster operation 
cannot be regarded as purely quasi - steady. We may postulate 
that the thruster operates in a quasi-steady mode during the high 
current portion of the current pulse, but departs from quasi- 
steady operation significantly during the fall and afterglow 
portions of the pulse. The net effect is to reduce the overall 
thruster efficiency for the pulse. 

An estimate of the ablated mass time history is given by the 
computational model for a 2 kJ shot with polyethylene in the 
quasi-steady mode (Figure 15). The model predicts a total ablated 
mass of 6.2 mg, in excellent agreement with the average measured 
value of 6.5 mg per pulse. At the end of the current puise, 
however, the model predicts that only 4.1 mg has been ejected. A 
significant fraction is thus exhausted at a much reduced velocity 
characteristic of the low temperature afterglow. 



COMPUTED CAPILLARY MASS LOSS 


ORIGINAL PAGE 19 
OF POOR QUALITY 
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VI. CONCLUSION 


A PET Thruster has been built and tested on a thrust stand, 
giving repeatable results over a series of shots. Exhaust velo- 
cities for polyethylene and Teflon propellants in the quasi- 
steady mode of operation are 27 km/sec and 20 km/sec, as measured 
by exhaust probes. These velocities are consistent with computa- 
tional predictions for these materials, assuming an adiabatic 
expansion to high Mach number In the conical nozzle. The perfor- 
mance data should be regarded cautiously, however, particularly 
with respects to the calculation of thruster efficiency. Future 
experiments are planned which provide for Impulsive calibration 
of the thrust stand, and for determination of the mass flow time 
history In order to provide a clear Interpretation of the 
ef f ici ency 


data. 



